Thermodynamics of the first— order vortex lattice 
melting transition in YBa2Cu307_5 
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Using the London approximation within the high field scaling regime, we calculate the jump in 
the specific heat Ac at the first-order melting transition of the vortex lattice in YBa2Cu307_4. This 
has recently been measured [A. Schilling et ai, Phys. Rev. Lett. 78, 4833 (1997)] and reported to 
be at least 100 times higher than expected from the fluctuations of field induced vortices alone. We 
demonstrate how the correct treatment of the temperature dependence of the model parameters, 
which are singular at the mean-field Bc2 line, leads to good agreement between the predictions of 
the London model and the size of the experimental jump. In addition, we consider the changes in 
the slopes of the magnetization A{dM/dT) and A{dAI/dH) at the transition. Using continuum 
anisotropic scaling theory we demonstrate the consistency of measurements at different angles of 
the magnetic field with respect to the crystal c-axis. 

PACS numbers; 74.60.Ec, 74.60. Ge 



I. INTRODUCTION 

In a recent paperEI we calculated the size of the jumps 
in the entropy and magnetization at the vortex-lattice 
melting transition within the London model. We have 
shown that, by correctly treating the temperature de- 
pendence of the model parameters, and knowing the vol- 
ume of the relevant fluctuation degrees of freedom, good 
agreement is obtained with the experimental results in 
bothYBa2Cn307_5 (YBCO) and BiaSrzCaiCuzOs (BiS- 
CCO) superconductors. The analysis contains only one 
unknown number, which can be taken from numerical 
simulations. The importance of these results is that they 
settle a controversy of recent years: it had been thought 
that j±|be observed jumps of order 1 ks per vortex per 
layeroa are incompatible with a simple melting scenario 
based on fluctuations of field induced vortices alone. Our 
analysis has demonstrated that there is no incompatibil- 
ity because the temperature dependence of the model 
parameters reflects the underlying microscopic degrees 
of freedom. This allows the correct size of jumps to be 
found without having to explicitly include extra fluctua- 
tions. 

Another characteristic of the melting transition is the 
jump in the specific heat capacity Ac — TA{ds/dT) = 
—T A{d^ g / dT^) , where g is the Gibbs free energy den- 
sity. Tphere have been careful measurements of this in 
YBCOpQ and in Ref. ^ it was claimed that the step in 
the specific heat is at least one hundred times too large 
to be explained by the extra fluctuations of the trans- 
lational degrees of freedom in the vortex liquid. The 
main purpose of this paper is to explain the size of the 
specific heat step using the London model, following the 
same approach of Ref. |l|. We emphasize that at the rel- 



evant melting fields in the YBCO system, the penetra- 
tion depth A(r) is much larger than the distance between 
vortices gq ~ {^q/BY^^ {B is the magnetic induction, or 
flux density, and $0 is the flux quantum) . In this regime 
the London model has simple scaling properties, which 
may be used to find the exact form of the jumps at the 
transition. 

In the next section we calculate the heat capac- 
ity within the London model, and compare our es- 
timated jump at the tra nsition with the experimen- 
tal values. In Sec. Ill we consider the jumps in 
the magnetization slopes A{dM/dT) = -A{d'^g/dHdT) 
and A[dM/dH) = -A{d'^g/dH^) and find results coi^ 
sistent with experimental values given by Welp et al.U 
We also take the opportunity to compare the jumps in 
the entropy and the specific heat between a system at 
constant external field (the experimental scenario) and 
a system at constant vorte x de nsity (the case for most 
simulations). Finally in Sec. |l^ we consider the effects of 
rotating the magnetic field away from the c-axis. p^is has 
been done in recent specific heat measurements. Qi 



II. THE HEAT CAPACITY IN THE LONDON 
MODEL 

We first consider an isotropic superconductor in the 
mixed state. Within the London approximation, the 
magnitude of the superconducting order parameter is 
taken to be constant, except within the vortex cores, 
which are assumed to be much smaller than the distance 
between vortices. These assumptions will hold at low 
enough fields below the upper critical field iJc2- The free 
energy of the system can then be expressed as aj-sum over 
pairwise interactions between vortex seginents,l3 i.e. 



J^L[{rp}] = ^51 / d-^i^-dr^ 



fJ.V 



-r„|/A 



(1) 



Disregarding the small distance cut-off, ^, there are two 
length scales in this problem: the London screening 
length A and the average distance between vortices oq. 
The energy scale per unit length Eq is equal to ($o/47rA)^. 
It is simple to include a uniaxial anisotropy into the 
model when the effective mass in the field direction is dif- 
ferent from the other two perpendicular directions: For 
a ratio of mab/rnc = £^, the lengths in the field direction 
are scaled by e. 

The experiments of interest are carried out at constant 
external field H . For the case of large field melting, as in 
YBCO, the magnetization is small and B k, H. In this 
regime we have ao <C A. and the London free energy takes 
a simple scaling formja 



^L[{r^}] ~ eao£o/L[{s^}], 



(2) 



where s^ = (x^, y^i, Zp/e)/ao are dimensionless posi- 
tion vectors and /l[{s^}] is a dimensionless functional 
independent of A. The thermodynamic properties of 
the system are determined by the partition function, 
Z = Trexp(— eao£o/L[{s^}]/r). Notice that one may 
think of T = T/eaoEo as an effective dimensionless tem- 
perature. From the standard thermodynamic relations, 
F = -T\nZ, S = -{dF/dT)B, and E = F + TS, the 
entropy is 



S 



T T \dT 



(3) 



The contribution S'o — {{T) — F)/T may be thought of as 
the configurational entropy of the coarse grained model, 
while the last term in (y) represents additional contribu- 
tions from the underlying microscopic degrees of freedom 
that appear in the temperature dependence of the model 
parameters. 

In the London model we must include the temperature 
dependence of the energy scale So{T) which we show in 
Ref. |l| to be important for an adequate determination of 
the jump in the entropy at the melting transition. Using 
the scaling of Eq. (H), the entropy jump isEl 



--"-£§)-- 1^ 



A^o- 



(4) 



The last form in Eq. M) uses the phenomenological tem- 
perature dependence of the penetration depth, A(r)^ = 
Ao/(l — P), where t = T/Tc is the reduced temperature. 
The length Aq is found from fitting this form to the slope 
of 1/A(T)^ as the zero field transition temperature Tc is 
approached,Q which for YBCO gives Aq =1300 A. 
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FIG. 1. The configurational average of the London free en- 
ergy for temperatures near melting, as measured in the simu- 
lations of Ref. 10 for a system of 81 lines, in the incompressible 
limit ao <C A. The dashed lines are straight line fits for the 
solid and liquid phases. 

The configurational entropy jump A5'o may be deter- 
mined from numerical simulai|ians. In the recent simula- 
tions of Nordborg and BlattertB the expectation value of 
/l[{s^}] is measured at different temperatures for a sys- 
tem at infinite A, the results of which we show in Fig. [^. 
This simulation uses the approximation of only including 
interactions within the same plane perpendicular to the 
external ifield, i.e., "retarded" interactions (in the Bose 
languagaa) are ignored. The volume of the system is 
V = LzNa^\fil2 with the number of vortices given by 
A^, and the system length in the field direction equal to 
Lz (periodic boundary conditions were used in all three 
directions). Inspection of Fig. n^ gives 



A5o 



A(/l> 



0.17iVi./ao, 



(5) 



where the transition occurs at r = t,„ — 0.088. Eq. (0) 
shows that the jump in configurational entropy is a 
constant per elementary degree of freedom: AS'o — 
0.17y/T4df, where Kdf = eao*&o/^ ~ eoo- The scal- 
ing factor (1 -|-t^)/(l — t^) in (0) then leads to a constant 
jump in entropy per vortex per superconducting layer. 
In Ref. d this was found to be AS'd « 0.4fcB when YBCO 
parameters were used, consistent with the experimental 
result measured in Ref. ||. 

The main point of this paper is to apply similar con- 
siderations to the heat capacity. This can be defined 
for changes at constant field Fl or at constant flux den- 
sity B. In the incompressible limit (which coincides with 
A/ao — > oo) these will be the same, and 



C, 



H 



Cf_ 



T 



dS_ 
dT 



dE 
df 



(6) 



(In the next section we will consider the difference 



Ch — Cb for a compressible system, and show that the 
difference between the jump in specific heat at constant 
B and at constant H is negligible on the YBCO melting 
line.) We again use the results of Ref. |l^ (shown here in 
Fig. nl) but this time to determine the jump in the heat 
capacity. 

Let us first assume that Eq is independent of temper- 
ature. Using (^), the heat capacity at constant density 
takes the form, 
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(7) 



Note that this "bare" heat capacity is simply a measure 
of the amplitude of energy fluctuations, which can be 
expected to jump at the transition from solid to liquid. 
It turns out that the first form is more accurate to use, 
once many temperatures have been sampled in the sim- 
ulations. The change in slopes in Fig. |l] gives the result 



AC. 



-BO 



OMNLjaQ. 



(8) 



This is interpreted as a jump in heat capacity of 0.38 ks 
per vortex degree of freedom (defined above). As 
Schilling et al. point out, this result is far too small to 
account for the jump in specific heat observed in their 
experiment. However, we have neglected the additional 
terms arising from the temperature dependent parame- 
ters, which we have seen give dominant contributions to 
the entropy jump close to Tc- For a general temperature 
dependence of the energy scale eo{T), we find that the 
heat capacity is. 



C, 



B 



T^ d'eo (/l) 
eo dT"^ T 



^-tMP- 



(9) 



Note that we need both results ^ and (||) to find the 
jump in Cb, although as Tc is approached and eo -^ 0, 
the last term in (0) dominates. 

In the following, we consider a more complex depen- 
dence for the energy scale of the vortex system cq (T) than 
was used in Ref. [^ and include a correction that accounts 
for the suppression in the superconducting density as the 
mean-field Bc2 line is approached. This i?c2 correction is 
important for a quantitative fit to London theory of the 
melting fine of YBCO at fields above 1-2 T. While good 
agreement for the entropy jump was found in Ref. be- 
tween the prediction of the uncorrected London theory 
and the experimental values, we will see that the spe- 
cific heat jump determined from (M is more sensitive to 
the exact form of the melting line. This is because the 
terms arising from the internal temperature dependence 
of £o(T) are more strongly diverging in this case than for 
the entropy jump in Eq. (^. 

In the limit B « Bc2, the system is described by 
Ginzburg-Landau (GL) theory, with an order p^u. 
rameter ?/; restricted to the lowest Landau level.tiJ 
The GL free energy takes the form Fql = 
/d3r[-a|V'P+/3|V'|* + 7|9V'/5z|2], where a{T,B) = 



ao(l - h) and h = 5/5^2 (T) (see for example Ref. H). 
Notice the dimensional reduction in the derivative term; 
there is a dijierging GL coherence length in the field di- 
rection only.Ej The only length scale remaining in the two 
perpendicular directions is the magnetic length ag. The 
order parameter scale is |'0oP — ck/2/3 and the conden- 
sation energy scale is 2a|'0oP = cP' / P- For this reason, 
the shear modulus haa_the behavior cge ~ (1 — 6)^ (as 
calculated by Labuschtfl). However, the tilt modulus at 
short wavelengths (equivalent to the supecfluid density in 
the 2-direction) has a linear dependencelij C44 ~ (1 — 5). 
(The correct thermodynamic limit of a constant C44 as 
kz —^ Q '\s not captured by the above LLL free energy, 
which is restricted to wavelengths below the penetration 
depth, k^ > 1/A. The full calculation in Ref. |l^ recovers 
this limit within GL theory.) 

Our intention is to use an effective pairwise interaction 
in (n[) that gives the London model at low fields, but with 
the correct elastic moduli in the high field limit. We take 
the scale of interactions to be £0 (2^) = ^o (T) ( 1 — 6) ^ . This 
then accounts for the suppression of the order parameter 
around both vortex segments in an interaction term. To 
account for the increasing stiffness in the field direction 
we scale the effective anisotropy e = e{\ — b)~^^^. Us- 
ing these effective parameters in the London model we 
find cgg oc 6(1 — 5)^, in agreement with the results of 
GL theory -at intermediate fields that have recently been 
calculated, t£l and C44 oc (1 — b). A similar extrapolation 
technique was used twenty years ago by Brandt .113 Con- 
verting the effective London system (with energy scale 
£0 and anisotropy £) to an isotropic one by scaling all 
lengths in the field direction, as in Eq. (g), gives the 
overall energy scale factor. 
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FIG. 2. London scaling fits to the melting line observed 
in Ref. 3 (circles) with (full line) and without (dashed line) 
Hc2 corrections. The Hc2 line we used in shown as the dotted 
line. The parameters used are given in the text. 



with £00 — (^o/47rAo)^. In the last line we have as- 
sumed the form BciiT) = Sc2(0)(l — i^) and written 
b = B/B,2{0). 

We now compare our results to the measurements of 
Schilling et alB^ and concentrate on the case where the 
field is directed perpendicular to the Cu02 layers. In 
the scaling regime, the melting line is described by a 
fixed value of t = Tm, i.e. Tm = gaoiB)eQ^{T„i, B). 
In Fig. we show our fit to the experimental melting 
line, where we find a value oi g ~ 0.21 (this corresponds 
to a Lindemann number of cl — 0.24 when the melt- 
ing line from the Lindemann criterion is written in the 
standard formcl of Tm ~ 2^y^^cj^eeoao)■ To obtain this 
fit we took the physical values: Tj, = 93 K, e = 1/8, 
dHc2/dT\T=T, = 1.8 TK-\ Ao = 1300 A. In Fig. | we 
also show our fit without including Bc2 corrections, which 
leads to a value of g ~ 0.13, or cl — 0.19. 

With ( |lO| ) the formula for the entropy jump in (Q) now 
becomes, 



t2[2(l _ i2)3 _ ^(1 _ ^2)2 _ ^2(1 + 2i2)] (/^) 



(l-i2_6)(l_i2) 



[ii-t^){i-t^-bW 



(11) 



A- 



The result for ACb in (gj) combined with (g) and (]5|) 
tells us the jump in the heat capacity per degree of free- 
dom (i.e., per volume T4df), which diverges as Tc is ap- 
proached. Indeed, we find that over the region of tem- 
peratures where measurements have been performed, the 
heat capacity jump is of order 100 ks per vortex degree 
of freedom. Because the melting field falls to zero, the 
volume of a degree of freedom also diverges in this limit, 
such that the heat capacity jump in a physical sample 
of fixed volume drops to zero, ACb oc (1 — t^), on ap- 
proaching Tc- In Ref. the results for the jump in heat 
capacity are expressed for a fixed volume. The specific 
heat is defined as the heat capacity in a mole of YBCO 
which fills a volume V^noi = 1.05 x 10~^m'^ (this is the 
volume of 6.02 x 10^^ unit cells). We therefore convert to 
the specific heat jump Acb = {Vinoi/V)ACB and make 
use of the fit to the melting line of Fig. || in Eq. (|g), 

ACbo - 0.381//14dt oc bI(\t). 



Note that the factor on the RHS now diverges at the 
Bc2 line rather than at Tc- In Fig. we compare our 
predictions for the entropy jump per vortex per layer, 
calculated with and without Bc2 corrections, to the ex- 
perimental values of Rcf. ^j. Interestingly, there is little 
difference between the two approaches in this tempera- 
ture regime. 




FIG. 3. The calculated entropy jump per vortex per layer 
(solid line). For comparison the results when Bc2 corrections 
are not included are also shown (dashed line) . Both curves are 
within the error bars of the experimental values from Ref. 3. 



Substituting ( fi(\ ) into (^ leads to 

[{1 - t^)il - t^ - b)]^ 




(12) 
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FIG. 4. The calculated jump in specific heat (solid line) 
with the units used in Ref. 4, and the measured points from 
that paper. In this case there is a significant difference from 
the calculation ignoring Bc2 corrections (dashed line) 

Our result for the specific heat jump per molar volume 
is shown in Fig. |4| along with the experimentally mea- 
sured points. This shows that we can explain the size of 
the specific heat jump within about 50% accuracy. We 
should not expect a much better agreement as we ex- 
trapolate the London approximation into a region where 
it will not be exact and also neglect the effects of dis- 
order. In fact, the measurements of the magnetization 
jump in Ref. [g show AB falling to zero at a temperature 
several Kelvin below Tc, and this is usually attributed 
to the effects of quenched disorder adaich may smooth 
the jumps at a first-order transition.li3 There may be a 
similar explanation behind the non-monotonic behavior 



in the experimental value of Ac. Nevertheless, the im- 
portant point here is that the size of the experimentally 
observed jumps in the specific heat are consistent with a 
model that only includes field induced vortices. 



III. THE MAGNETIZATION AND ITS 
DERIVATIVES 

The aim of this section is to use our results for the en- 
tropy and specific heat jumps from the previous section 
to calculate the jumps in the derivatives of B (or in the 
magnetization M = {B — H)/A-k) with respect to H and 
T. We first introduce the Clausius-Clapeyron relation, 
and then describe the conditions under which the jump 
in entropy in a system at constant external field H is the 
same as the entropy jump in a system of fixed flux density 
B. We need these conditions to hold in order to justify 
our application of the results of Section y, which we cal- 
culated at constant B, to the experimental results, which 
are measured at constant H. After ensuring that these 
conditions are fulfilled for the melting line of YBCO, we 
then calculate the jumps in the thermal compressibility 
A{dB/dT)H and in the susceptibility A{dB/dH)T. Fi- 
nally we will be in a position to describe the conditions 
for the specific heat jump to be approximately the same 
in a constant B system compared to that in a constant 
H system, and then verify that these conditions are also 
satisfied. 

At the first-order transition of a system at constant 
field H, the continuity of the Gibbs free energy G — 
F — {BH/A-k)V leads to the Clausius-Clapeyron equa- 
tion 



As = 



1 dH„ 



AB, 



(13) 



relating the jump in flux density AB to the melting 
line Hm(T) and the jump in the entropy density As. 
In the previous section we calculated the jumps in en- 
tropy and heat capacity for a transition at fixed B, us- 
ing the scaling form of the London model [see Eq. (0)]. 
An important difference between the melting transition 
at constant B and at constant H is that in the former, 
coexistence occurs over a finite temperature region, see 
Fig. g. The entropy jumps at constant H and constant 
B will not be the same: We must take into account the 
different values of B on either side of the constant H 
transition, AB\h = B^{H,Tm) - B''{H,Tm), and also 
the width in temperature of the constant B transition, 
AT\b = Tl^iB) - T^,-XB) [the symbol s denotes the sohd 
(lattice) phase, while the symbol f is for the fluid (liquid) 
phase]. The entropy jump at constant B is defined as 
s^iB,Tl^) - s''iB,T^) as shown in Fig. |, which to first 
order in AT is 
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FIG. 5. Schematic phase diagrams in the H-T plane, 
showing the phase boundary Hm{T), and the B-T plane, 
showing a region of coexistence of the fluid and lattice phases 
between the two boundaries B'^{T) and B (T). a) A path 
is shown for increasing temperature at constant external H, 
crossing a sharp phase transition at Tm{H). In the B-T plane, 
the coexistence region is crossed at constant T, with a jump 
in induction AB\h- b) The path for increasing temperature 
at constant B is shown. While traversing the coexistence re- 
gion, the path remains on the phase boundary Hm{T) in the 
H-T plane. There is a well defined width of the transition, 

Arjfl. 




FIG. 6. Schematic graph of the entropy as the melting 
transition is crossed at constant flux density B. The en- 
tropy jump is defined across the whole width of the transition, 

As\b = s\B,tI,)-s%B,T^). 



As\B^s'iB,T:,, 



^{B,t:,,) + ^AT\b. 



(14) 



We can now relate the entropy jump at fixed H to that 
at fixed B, 



As\H^s'iH^,T)-s'{H^,T) 



As|s= + ^ 



1 fdH^ 



An \ OB J ^ \ dT y jj 



'dB 
dB^ 



T 

ASIh 



H 



-f^T\B, 



(15) 



where in the last hne we have used the MaxweU relation 



ds_ 
'dB 



4^ ['df 



4^ ['dB 



dB 
'df 



H 



(16) 



Equation ( p^ ) shows that there are two terms that lead 
to corrections in the entropy jump at fixed B compared 
to fixed H. The physical explanation is that the latent 
heat at fixed H includes the magnetic work done to in- 
crease the flux density by AB\h, while at fixed B there 
is an extra heat increase due to the temperature increase 
across the transition AT\b- Our results for the entropy 
jump at constant B will only be valid at constant H, 



A 



s\h 



A.s|„= 



B=B{H) 



i_ /dip' 

Att V dB 



if the two conditions 

^AT|b^„«A.|^ (17) 

J^^) AB\h<^As\h (18) 



are fulfilled. 

To apply the results of the previous section, we must 
check the above conditions for the case of vortex-lattice 
melting in YBCO. We first consider the term arising from 
the temperature increase across the constant B melting 
transition. By writing AB\h ~ {dBm/dT)AT\B (see 
Fig. H), and using the Clausius-Clapeyron relation (13), 
we can express the condition (n^) in the form, 



c4 1_ fdBrn^ 

T Att\ dT 



(19) 



Using the numerical results presented in the previous sec- 
tion, the total heat capacity at the YBCO melting tran- 
sition is of order cb/T sa 10 Jm~'^K~^. From Fig. || 
we have dHm/dT w — 0.7TK~^ on the melting line of 
YBCO, which gives (dHm/dT)^ /An « 4 x 10^ Jm-^K'^. 
Therefore this condition is well satisfied on the YBCO 
melting line. Finally, we consider the term involving 
the jump in flux density at the transition with fixed H. 
If we insert the Clausius-Clapeyron relation ( [l^ ) into 
the RHS of condition (p^ ) we find that the condition 
is equivalent to [dB^ /dT)H ^ dHm/dT, which is to 
say that the thermal compressibility of the vortex sys- 
tem {dB^/dT)H = An{dM^/dT)H must be smaller than 
the slope of the melting line. This is satisfied as long 
as the magnetization (which is of order the lower criti- 
cal field Hci) is less than the melting field, Hd ^ i?,„. 



which is the case for vortex lattice melting in YBCO. We 
have therefore confirmed in our regime of interest (the 
large field scaling limit) that the same entropy jump is 
found for a transition at constant H and at constant B. 
If we also take Bm{T) — Hm{T), we can calculate AB in 
the scaling regime using the Clausius-Clapeyron relation 
(O). This reasoning was used in Ref. ^ and excellent 
agreement was found between the London model predic- 
tionS|^d the measured values of AB in YBCO by Welp 
et aim 

We now calculate the jumps in the derivatives of the 
magnetization. Consider a quantity X{T,H) which is 
discontinuous at the melting field Hm{T) with a jump 
AX{T). The total temperature derivative of the jump in 
X (defined along the melting line) is then 



dAX 
dT 



dX 
'df 



H 



dHm.fdX 
dT [dH 



(20) 



Note that inserting X ^ G, the Gibbs free energy which 
is continuous, simply gives the Clausius-Clapeyron equa- 
tion. Applying ( po| ) to the entropy jump by taking 
X = s, and using the Maxwell relation An{ds/dH)T = 
{dB/dT)H, we find 



1 dHm 

4^ dT 



dB\ 
dT 



AcH dAs 



H 



T 



dT 



(21) 



This allows us to find A{dB/dT)H once we know the cor- 
responding discontinuities in the calorimetric quantities 
ch and s. In the SQUID measurements of the magne- 
tization in Ref. o, an experimental value is found at a 
field oiH = 4.2T of A{dB/dT)H = AnA{dM/dT)H ~ 
0.2 GK~^ (the accuracy is only of order 50%). We con- 
vert from magnetic energy units to calorimetric units 
using lergcm"^ = O.IJ m^'^ = 0.011 mJmoP^ for 
YBCO. For the specific heat jump we take Ach ~ 
Acb = 1 rnJ moP^K"^ from Fig. ^. We calculate 
As = AS/V from Eqs. (||) and (|) giving dAs/dT = 
-0.35mJmor^K-2 at a field of Hm = 4.2T. Insert- 
ing these results into (|2^) along with the melting slope 
of dHm/dT = -0.7TK-1, we find A[dB/dT)H = 
0.23 GK^^ consistent with the experimental value at this 
field. 

We now use two further relations between these 
jumps which will allow the jump in the susceptibility, 
A{dB/dH)T, to be calculated. The total derivative with 
respect to temperature of ( |l3| ) may be written as 
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This equation was used in Ref. ^ to successfully test the 
thermodynamic consistency of the observed specific heat 
jump with the magnetization measurements of Welp et 



al.a Combining this with (|21|) gives an equation for the 
jump in {dB/dH)T in terms of calorimetric quantities 
only, 
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Note that this may also be derived by inserting X = 
B in Eq. (^ and using the Clausius-Clapeyron re- 
lation. We find this jump along the melting line of 
YBCO to be of order A{dB/dH)T « lO-^G/Oe, con- 
sistent with the assumption that B k, H. However, 
we must be careful about ignoring this term. For in- 
stance, in the scaling regime, and from experimental val- 
ues in YBCO, the temperature derivative d{AB)/dT — 
A{dB/dT)H + {dHm/dT)A{dB/dH)T is negative but 
the jump A{dB/dT)H is positive and therefore smaller 
in magnitude than {dHm/ dT)A{dB / dH)T- 

We now determine whether the jump in the heat ca- 
pacity at constant iJ, as measured in experiments on 
YBCO, may be approximated by the jump in heat ca- 
pacity at constant _B, which was calculated in Section ||. 
As in our earlier consideration of the entropy jump, we 
must include the change in B at a constant H transi- 
tion as well as the finite temperature width of a con- 
stant B transition. However, there is an extra differ- 
ence in that, while the entropy is a function of state, 
S{H,T) = S[B{H,T),T], the heat capacity is defined 
differently for changes at constant B or at constant H, 
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where the last line again uses the Maxwell relation ( |16[ ). 
With the same reasoningthat lead to (|lj) and using the 
notation defined as Fig. ^ we have, 



the total Ch ~ cb is rigorously positive for thermody- 
namic stability, the change in this quantity may be of 
either sign. Combining this equation (p&) with (pS) gives 
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We can group the corrections to Acb\b mto two con- 
tributions. The first correction term in ( p7| ) can be 
ignored compared to the total specific heat jump if 
{dc^jj / dT) AT\b ^ Acs. After some manipulations with 
the Clausius Clapeyron relation it may be written as 
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(28) 



which we find to be satisfied in the scaling regime from 
the simulation results. A second condition is found from 
the three remaining correction terms in (Ujh after drop- 
ping the term oc AT\b- Setting {dH/dB)^ — 1 and using 
the identity {dcH/dB)T = -T{dH/dB)Tid^B/dT^)H 
[which also arises from the Maxwell relation dl9)], gives 
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A direct comparison of this equation with the differen- 
tiated Clausius-Clapeyron relation (Ea) shows that we 
can take Ach\h ~ Acb|b as long as \dH,n/dT\ » 
{dB/dT)H, which is the same condition we found for 
comparing the entropy jump at constant H versus con- 
stant B, and which is valid as long as Hm ^ Hci- There- 
fore, our calculation of Acb at constant B is justified in 
the application to vortex lattice melting in YBCO. 
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We can also find the jump in the difference [ch — cb) 
from (bj) by expanding to first order in the jumps. 
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Our previous results show that the first term on the RHS 
is negative, while the second term is positive. Although 



IV. DIFFERENT FIELD ANGLES 

In Ref. the heat capacity in YBCO was measured 
both with the magnetic field parallel to the c-axis, H || c, 
and parallel to the superconducting layers, H _L c. More 
recent experiments have been performed for a range of 
angles between these two limits.Ll In the scaling regime, 
the effects of rotating the field may be simply understood 
using the anisotippic scaling rules of Blatter, Geshken- 
bein, and Larkin.ll3 The basis of these scaling rules is that 
the anisotropic system can be transformed to an isotropic 
one by scaling the lengths along the c-axis. For a mag- 
netic field enclosing an angle 9 with the c-axis, a given 
quantity Q then follows the scaling law 



Q{9, H, T, Xab, e) = sgQiseH, T/e, A), 



(30) 



where Q is the corresponding quantity in the isotropic su- 
perconductor (with A = Xab), and e^ 



-^ sin^ ' 



The factor sq depends on the quantity Q. For the scalar 
quantities volume and energy it is independent of the 
angle, sq — e, but for magnetic fields it is equal to 
SB = l/ee- 



An initial consequence of ( |30[) is that, at a fixed tem- 
perature, the magnetic field at which the system has the 
same physics scales as Hx{0) = Hx/se- This implies 
that the melting line is raised to higher fields as the an- 
gle 9 is increased. 



H„-,iT,9)^HUT,0)/sg. 



(31) 



The same scaling form applies to the upper critical field 
Hc2{T, 9). These fields are lower for H || c than for H _L c 
by a factor s. Furthermore, if we take the quantity Q in 
( po| ) to be the energy, entropy or heat capacity of a fixed 
volume (e.g. the volume of the crystal) then after scaling 
the magnetic field as in (pW this quantity will be inde- 
pendent of the angle 9. Therefore, if we write the entropy 
jump as a function of melting temperature, As(T, 9), the 
following relation holds. 



As[Tr,^iH,9),9] = As[TmiseH,0),0]. 



(32) 



We can also replace s by c in this equation. Therefore the 
measured jumps in the specific heat and in the entropy 
per molar volume, when plotted against the melting tem- 
perature, will be identical for the different angles. This 
prediction is fulfilled to within experimental error in re- 
cent measurements by SchillingJII 

To understand this result in more detail, we now 
consider the results for the jumps (Q) and (p^). No- 
tice that both of these results are of the form AX — 
f{T/Tc,B/Bc2)V/Vcdi. The scaling rule ^ tells us 
that at fixed T, the reduced variables t = T/Tc and 
b = Bjn/Bc2 at melting are independent of 9. We there- 
fore concentrate on the angular dependence of V^df- In 
the isotropic system, V^df = Og. However, the three per- 
pendicular dimensions of V^di will scale differently with 
9. If the field lies in the y-z plane then r^^ = I =|-p.Oi 
whereas the lengths in the plane of rotation scale as,l 



(33) 



Therefore, at fixed temperature, but scaling the magnetic 
field as in Eq. (^), we find that the volume per degree 
of freedom stays constant, 
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So we see that it is because the volume per vortex de- 
gree of freedom is independent of 9 at the melting field 
of a fixed temperature that As{T) and Ac(T) are also 
independent of 9. 



We now consider the case of fixing the magnetic field, 
and comparing the jumps at different melting tempera- 
tures as the angle of the field is rotated. In this case, the 
field scales to a different value, H — sgH, in the isotropic 
system for different values of ^. The melting temperature 
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(35) 



This equation differs from (pVj in that we need to know 
the full melting curve in the isotropic case to find the 
angular dependence of the melting temperature. Sim- 
ilarly, to determine the angular dependence of the en- 
tropy jump as a function of H from (B3), we need to 
know the full function As(T, 9 = 0) and the melting 
curve Tm{H, 9 = 0). However, in the case of no Bc2 cor- 
rections we have an analytic form for the melting curve, 
Bm{T,9) ex (1 — t'^)'^/ee, and in this limit we can find 
the angular scaling at the same values of H. When we 
stay on the melting line at the same Hm, the volume 
Vcdf{9) = eV'cdf scales as 
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dt 
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3/2 



(36) 



If the jump in the entropy and the heat capacity were 
constant per vortex degree of freedom then this equa- 
tion implies that the jumps would scale proportional to 
Eg' , as is stated in Ref. ^. However, we showed in Sec- 
tion H that the temperature dependence of eo{T) gives 
a non-trivial temperature dependence to the jumps per 
vortex degree of freedom. In the limiting case of no 
Bc2 corrections, the jump in entropy density is As ex 
(1/K!df)(l — t^)^^ ex BmSe and the jump in heat capac- 
ity is Ac ex (1/V;df)(l - i^)"^ oc (S,„ee)i/2. This leads 
to the angular scaling 



As{H,9) = As{H,0)ee 
Ac{H,9) =Ac{H,0)ey^. 



(37) 



For the general case when Bc2 corrections are included, 
there is no such simple form. However, the more gen- 
eral angular scaling of the jumps ( p^ ) still holds and it 
is straightforward to find As{H, 9) and Ac{H, 9) numer- 
ically. 



V. CONCLUSIONS 

We have not in this paper considered the case of the 
strongly layered high-Tc superconductors such as BiS- 
CCO, where the low melting field Bm < ^o/^^ means 
that the scaling form of (||) does not apply. Still, in Ref. [|, 
we successfully determined the magnetization jump in 
BiSCCO using dimensional estimates, together with a 
Lindcmann analysis of the melting line. The reason this 
works, while simple estimates for the entropy jump do 
not, is that the induction is a derivative of the Gibbs free 



energy with respect to H (rather than T as for the en- 
tropy) and the H dependence of the vortex parameters 
is much weaker than the T dependence. The entropy 
jump for BiSCCO could then be found by combining our 
result for AB with the melting line and the Clausius- 
Clapeyron relation ([l3|). It is not possible to simply fol- 
low this procedure to find Ac using the relation (|2^), as 
the simple estimates will be incorrect for A{dB/dT)H- 
It is also interesting to consider whether the conditions 
for the calorimetric jumps to be the same in a constant B 
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system and in a constant H system, derived in Sec. [II 



hold at these low fields. The melting field in BiSCCO is 
of the order of the bulk lower critical field Hd, which is 
also the size of the magnetization: Hm ~ Hd ^ AttM. 
Therefore the slopes dH^/dT and {dBjdT)H wiU be of 
a similar size so that the condition (nS) is not fulfilled 
and the RHS of ( |29| ) is not small. In experimental mea- 
surements on BiSCCO (at constant H) there will be a 
significant contribution to the entropy and specific heat 
jumps that is not present in a system of fixed flux density 
B. 

To summarize, we have calculated the specific heat 
jump at the melting transition of the vortex lattice for 
YBCO, using the scaling form of the London model at 
high fields and numerical results from recent simulations. 
We find a jump with the same magnitude as the measured 
experimental values. We have used analogous relations 
to the Clausius-Clapeyron equation to find the jumps in 
the magnetization slopes, and these are also consistent 
with measured values for YBCO. A careful analysis was 
made to make sure that the jumps in the entropy and in 
the specific heat are the same for a system at constant 
B and at constant H . We have constructed the relevant 
conditions and have verified that they are satisfied at the 
vortex-lattice melting line in YBCO. We have analyzed 
the changes to the jumps when the field angle is rotated 
with respect to the crystal and thereby explain the re- 
cently measured angular scaling of the entropy and spe- 
cific heat jumps in YBCOlO In conclusion, we find that 
the London model with temperature and field dependent 
parameters gives a consistent-picture of the first-order 
transition observed in YBCOBj 
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